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ABSTRACT 

Spectroscopic observations of Ho and H/3 emission lines of 128 star-forming galaxies in the redshift 
range 0.75 < z < 1.5 are presented. These data were taken with slitless spectroscopy using the G102 
and G141 grisms of the Wide-Field-Camera 3 (WFC3) on board the Hubble Space Telescope as part 
of the WFC3 Infrared Spectroscopic Parallel (WISP) survey. Interstellar dust extinction is measured 
from stacked spectra that cover the Balmer decrement (Ha/H/3). We present dust extinction as a 
function of Ha luminosity (down to 3 x 10 41 erg s _1 ), galaxy stellar mass (reaching 4 x 10 8 M ), 
and rest-frame Ha equivalent width. The faintest galaxies are two times fainter in Ha luminosity 
than galaxies previously studied at z ~ 1.5. An evolution is observed where galaxies of the same Ha 
luminosity have lower extinction at higher redshifts, whereas no evolution is found within our error 
bars with stellar mass. The lower Ha luminosity galaxies in our sample are found to be consistent 
with no dust extinction. We find an anti-correlation of the [O III]A5007/Ha flux ratio as a function 
of luminosity where galaxies with Ln a < 5 x 10 41 erg s _1 are brighter in [O III]A5007 than Ha. This 
trend is evident even after extinction correction, suggesting that the increased [O III]A5007/Ha ratio 
in low luminosity galaxies is likely due to lower metallicity and/or higher ionization parameters. 
Subject headings: dust, extinction — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Star-forming galaxies are characterized by emission 
lines from gas heated by young stars in H II regions. The 
luminosity of these nebular emission lines such as Ha or 
[O III] provides reliable information concerning galaxy 
star-formation rates (SFRs) and gas-phase metallicities. 
Yet the dust present in the interstellar medium strongly 
attenuates rest-frame ultraviolet and optical fluxes in 
a wavelength-dependent manner. Therefore, inferring 
physical properties (e.g., SFRs) from emission- line lumi- 
nosities requires accurate dust corrections. 

In practice, the most reliable technique to estimate in- 
terstellar extinction is to measure the flux ratio of two 
nebular Balmer emission lines such as Ha/H/3 (i.e., the 
Balmer decrement). Since the value of the Balmer decre- 
ment is set by quantum physics, any deviation from this 
expected value may be attributed to dust extinction (for 
a fixed electron temperature). Other emission-line ra- 
tios, such as H/3/H<5 are useful as well, but are generally 
more difficult to measure because the higher order lines 



are weaker and more affected by stellar absorption. 

The determination of dust extinction from the Balmer 
decrement has been shown to be a very successful tech- 
ni que in the local Un iverse since the first statistical work 
by Kennicutt (1992 1. These results have been improved 
upon by the large amount of optical spectra provided 
by the Sloan Digital Sky S urvey (SPSS), which wer e 
analyzed in this context by [Brinchman n e t aL"| ( 2004 ) 
Moustakas et al.| ( |2006 [), and|Garn & Best| (|2010p^ow 



ever, the simultaneous detection ot Ha and tip in higher 
redshift galaxies, which are shifted to the near-infrared 
(IR), is a difficult task with current instrumentation on 
ground-based telescopes. Indeed, Balmer decrements in 
high-redshift galaxies have been spectrosc opically mea- 



sured only in a few indivi dual cases (e.g., Teplitz et al. 



2000 |Hainline et al.||2009[ ). 

Other recent, although less robust, 
Balmer decrement are available at 
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estimates of the 
z ~ 0.5 for larger 
galaxy samples b ased on a comb ination of spectroscopy 
and photometry ( |Ly et aL| [2012). There are other dust- 
extinction estimates at z ~ 1.47 from narro w-band pho- 
tome try of the Ha and [O II] emission lines 
but these measurements are affectec 



2012) 



Sobral et al. 
by the un- 
known metallicities of the galaxies. At high redshifts 
(z > 2) the ultraviolet (UV) s lope is typically used to de- 
termine dust extinction (e.g.,|Stanway et al.|2005| Hathi 
et al.|[2008l |Reddy et al.||2010| |Finkelstein et al'||2011[ 



Wilki ns~et al.||20lT " Bouwens et al.||2012| [Hathi el al.| 

2012), where access to the Balmer lines is difficult. The 



relationship between UV slope and dust reddening de- 
pends both on the dust properties and distribution of 
the dust with respect to the stars. While the local cor- 
relation between UV slope and dust attenuation appears 
to hold for typical star-forming galaxies up to z ~ 2 
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(e.g., |Reddy et aT1|2010[ |2012[ ), the Balmer decrement 
offers a complementary method of deducing dust atten- 
uation, and can illuminate any differ ences in the nebular 



versu s stellar extinction in galaxies ( |Calzetti et al.||1994 
20001 



Previous works on dust extinction have focused on 
a rel atively narrow dyn amical range of Ha luminosity 
(e.g., Sobral et al.|2012" ). Furthermore, these works have 
been generally sensitive to the brightest star-forming 
galaxies, which tend to have larger reddening than fainter 
galaxies. As a result, the common procedure to cor- 
rect for dust extinction in high-redshift galaxies is to as- 
sume an extinction of one magnitude at the wavelength 
of Ha for all galaxy luminosities. This assump tion is 
based on t he mean extinction of local galaxies ( |Kenni-| 
cutt 1992). However, this local average has a signifi- 
cant scatter, which may potentially affect results derived 
from the dust correction. An additional issue is that dust 
extinction as a function of bolometric luminosity is ob- 
served to evolve with redshift such that galaxies of the 
same luminosity are more extinguished at lower redshift, 
due mainly t o the metal enrichment of the interstellar 
medium (e.g., IReddy et al.||2006HBuat et al.||2007||Bur- 



garella et al.||20U7 


Bouwens et al.||2 


M\ IReddy et al. 


2010 Buat et al.||20: 
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2012]). 



Telescope (HST) has now made it possible to carry out a 
statistically significant analysis of the Balmer decrement 
in high-redshift galaxies. The grism spectroscopy mode 
of the WFC3 instrument offers several advantages over 
ground-based spectroscopy. First, grism spectroscopy is 
slitless, so there are no concerns of variable slit losses 
at Ha and H/3 due to varying seeing, alignment er- 
rors, and atmospheric refraction. Second, the WFC3 
grisms provide a continuous wavelength coverage over 
0.85 < A < 1.65 fim, which is unaffected by hydroxyl 
(OH) airglow emission and the high opacity of the atmo- 
sphere at certain wavelengths (i.e., line measurements 
are obtained regardless of whether they lie within atmo- 
spheric windows or between sky lines). Third, the flux 
calibration with WFC3 is more accurate than ground- 
based instruments, as it is unaffected by variable telluric 
absorption from Earth's atmosphere. Lastly, the WFC3 
grisms obtain unbiased samples, whereas ground-based 
spectroscopy typically requires a preselection in order to 
assign slits. Furthermore, longslit followup at higher red- 
shifts has mostly focused on brighter galaxies. 

In this paper, interstellar dust-extinction of star- 
forming galaxies is statistically investigated in the red- 
shift range 0.75 < z < 1.5 over two orders of magnitude 
in Ha luminosities. Balmer decrements are measured us- 
ing spectroscopy and analyzed down to an Ha flux limit 
of approximately 5 x 10~ 17 erg s _1 cm~ 2 . These are 
galaxies a factor of approximately two fainter than ga lax- 
ies previously studied at z ~ 1.5 ( Sobral et_al.|2012 |. 

The paper is organized as follows. In *|5J we present 
our galaxy sample and the emission-line extraction tech- 
nique. Section [3] explains the necessary framework of 
dust extinction used in this analysis. The results are 
shown in Q Section [5] provides a discussion of our re- 
sults, their robustness, and a comparison with previous 
works. Finally, we summarize the main findings from our 
analysis in 

Throughout this paper, a standard cold dark matter 



cosmology is assumed, with matter density Q m = 0.3, 
vacuum energy density tt\ — 0.7 and Hubble constant 
H = 70 km s" 1 Mpc" 1 . 

2. GALAXY DATA SET 

2.1. Observations and data reduction 

The data set used in this work is part of the WFC3 In- 
frared Spectroscopic Parallel (WISP) survey. The WISP 
survey offers the opportunity to study the spectral prop- 
erties of galaxies in a way that is unbiased by target se- 
lection. It also has unprecedented sensitivity to galaxies 
on the faint-end slope of the galaxy luminosity function 
at redshifts where th e peak of the star f ormation in the 
Universe occurs (e.g. 



Reddy et al. 2008) 



A detailed description of the WISP progr am includin g 



technical details can be found in Atek et al. (2010 2011 ) 
Briefly, the WISP program is obtaining slitless and low 
resolution near-IR grism spectroscopy of more than 200 
independent and high-latitude fields by observing in the 
pure parallel mode of the HSTj^J The data are taken with 
the infrared channel of the WFC3 in a total of ~ 730 
orbits using two grisms: G102 (0.80 < A < 1.17 /im, R 
~ 210) and G141 (1.11 < A < 1.67 /im, R - 130). 

The WISP data also include direct-imaging photom- 
etry with WFC3/IR and UVIS channels in the F475X, 
F600LP, F110W and F160W filters. IRAC 3.6 yum pho- 
tometry from the Spitzer Space Telescope is available for 
a s ubse t (~ 47%) of the initial galaxy sample described 
in j pT2l 

The optical and near-IR photometry is obtained in 
matched isophotal apertures defined by the F110W imag- 
ing. The galaxies are first detected with SExtractor 



(Bertin & Arnouts 
ing the position in 



1996) in the F110W band, comput- 
ihe sky and the aperture parameters 
(i.e., semi-major axis and positional angle of the ellipse). 
These parameters are used for our own routine to ex- 
tract the fluxes as follows. All four HST bands are ex- 
tracted in the same elliptical apertures defined in the 
F110W images. The local background for each aperture 
is calculated within a square annulus using a mean with 
3<r clipping to determine the representative background 
per pixel. Then, the total background-subtracted flux 
is computed within the aperture by adding all the aper- 
ture pixels. The flux uncertainties are computed using 
Monte Carlo simulations. First, for each aperture pixel 
we randomly vary the flux according to a normal distri- 
bution with a width provided by an error map. Second, 
we repeat this process 1000 times, extracting the ran- 
domized flux in each realization. Finally, from the sim- 
ulated flux distribution, the ler uncertainty is calculated 
as the width of a Gaussian fit. Further details on the 
photometric analysis will be found in Bedregal et al. (in 
preparation) . 

The Spitzer IRAC 3.6 /im photometry is extracted in 
fixed apertures of 2.8 arcsec in radius. Then, this pho- 
tometry is corrected to total fluxes assuming that the 
galaxies are point source^ 10 ! This is a fair assumption as 
the Spitzer point-spread function is significantly larger 
than the typical WISP aperture defined in F110W. The 

9 In the pure parallel mode one or more instruments sample 
the HST focal plane while the prime program observes its planned 
target according to its desired visit schedule. 

10 http://irsa.ipac.caltech.edu/data/SPITZER/Enhanced/Imaging/ 
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Figure 1. The 5<r line sensitivity shown in |Atek et al.| l |2010| l 
is recast here as H/3 luminosity (dashed lines) and ria luminosity 
(solid line) as a function of redshift. This sensitivity corresponds 
to a line covering an area of 15 pixels. The blue-dashed line corre- 
sponds to the H/3 line when it is detected by the G102 grism and 
the red-dashed line when H/J is detected by the G141 grism. The 
integration time is 6400 s for the G102 grism and 2700 s for the 
G141 grism. These exposure times are typical for the observations 
used in this study. The discontinuity in the H/3 luminosity reflects 
the fact that H/3 falls either in the G10 2 gri sm or in the G141 grism 
depending of the galaxy redshift. See §2.2| for details. 

relatively small IRAC aperture was chosen in order to 
minimize any strong contamination from brighter neigh- 
bors. 

The direct imaging is corrected for Galactic dust ex- 
tinction given the color excess from Schlegel et al. fl!998|) 
and a ssuming the extinction curve ot |(Jardelli et al.| 
([1989). 



2.2. Galaxy sample 

Our initial sample contains 312 star-forming galaxies 
in 17 different fields in the redshift range 0.75 < z < 1.5, 
where both Ha and H/3 fall simultaneously in the WISP 
spectral coverage. In the overall spectrum, our line finder 
flags potential emission lines as significant peaks above 
the continuum. The emission lines are identified and val- 
idated by eye by two astronomers, and then the redshift 
is assigned based on the centroid of the most significant 
line, which is usually Ha (Colbert et al., in preparation). 
We include in the final sample only galaxies with no ex- 
ternal contamination (i.e., galaxies with no overlapping 
spectra). This condition removes 49% of the galaxies in 
our initial sample. This method assures that we accu- 
rately measure both continuum and emission-line fluxes 
that are not affected by contamination. We also make 
sure that only galaxies with secure redshifts are included 
in the final sample (that is galaxies where it is clear that 
the selected emission line in the spectrum is Ha and 
not [O III]). This contamination is avoided by looking at 
other emission lines at their expected wavelength. This 
extra condition removes another 5% of galaxies in our 
initial sample. Our final sample contains 128 galaxies 
(with 60% of these galaxies detect ed by IRAC 3.6 tx m). 

The WISP 5tr line sensitivity (|Ate ket al1|2010[ ) re- 
cast in luminosity is shown in Figure [T] as a function of 
redshift. This sensitivity corresponds to a line covering 
an area of 15 pixels on the detector. The integration 
time shown in Figure [l] is 6400 s for the G102 grism and 
2700 s for the G141 grism. These exposure times are typ- 



ical for the observations used in this study. The red-solid 
line is the line sensitivity for Ha, which for our redshift 
range is always detected in the G141 grism. However, 
H/3 is observed either by the G102 grism (for lower red- 
shift galaxies) or by the G141 grism (for higher redshift 
galaxies). These H/3-line sensitivities are shown in Fig- 
ure [T] with blue-dashed line for the G102 grism and with 
red-dashed line for the G141 grism. The H/3 line appears 
in both grisms at z ~ 1.3. At this redshift, and given the 
typical relative exposure times in each grism, the G141 
is approximately a factor of 2.5 less sensitive than G102 
in detecting H/3. 

2.3. Emission-line extraction 

A non-parametric method is developed to subtract the 
continuum independently from each grism. Because we 
know the exact wavelength of the emission line (from the 
automatic line detection code), we interpolate over the 
lines (5 pixels on either side). Then, a median filter is 
run with a 40 pixel window over the spectrum. This is 
the model of the continuum. These model continua are 
subtracted from the spectra. 

Each spectrum is later divided into different regions 
in order to extract the emission-line parameters such as 
fluxes and equivalent widths: a bluer part containing H/3 
and the [O III] doublet, and a redder part containing Ha 
and the [S II] doublet. In the bluer region, a function 
consisting of three Gaussians (one for every line in the 
region) plus a straight line (modeling the small residual 
continuum) is fit using a \ 2 minimization. In the redder 
region, since our resolution is not high enough to resolve 
the two [S II] emission lines, only two Gaussians (one for 
Ha and another one for the [S II] doublet) plus a straight 
line are fit. The [N II] doublet cannot be resolved from 
Ha at our resolution. The [N I I] contribution to the Ha 
fluxes will be discussed in §2.6| We impose multiple con- 
straints to the line fitting. Specifically, the full width half 
maximum of the Gaussians is constrained to a maximum 
expected size of 10 pixels (~ 1.3 arcsec). The wavelength 
solution has an uncertainty of ±30 A. Therefore, in the 
fitting procedure the centers of the emission lines are al- 
lowed to vary within that wavelength range. The ratio 
between the fluxe s of the two [O III] lines is set to 3.2 



(Osterbrock 1989). We note that in order to avoid ob- 
taining local minima in our fits, we use a grid of initial 
conditions to find the global minimum. 

Finally, Balmer decrements are computed as the ratio 
of the Ha and H/3 line fluxes. The uncertainties in the 
Balmer decrements are then calculated by propagating 
the uncertainties of each of the individual line fluxes, 
which are obtained from the \ 2 fits. 

2.4. Stacking of the galaxy spectra 

The individual measurements are not always good 
enough to detect the Balmer decrements. The Ha emis- 
sion line is detected in all galaxies by construction, while 
H/3 is detected at more than 3c in 26 galaxies (this is 20% 
of the final sample). Therefore, to improve the statistics, 
we next stack all the continuum-subtracted spectra. To 
stack these galaxy spectra, we first shift each spectrum to 
the rest-frame, interpolating the values to a finer, com- 
mon wavelength grid. Different stacked spectra are built 
selecting the galaxies according to their Ha luminosity, 
stellar mass, and Ha equivalent width (see Figure [2]). 
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Figure 2. Composite spectra in the redshift range 0.75 < z < 1.5. The top panels show the stacked spectra for three bins in Ha luminosity, 
the middle panels show the stacked spectra for three bins in stellar mass, and the bottom a panels show the stacked spectra for three bins 
of rest-frame Ho equivalent width. The number of stacked galaxies is 43, 43, and 42 in the low, middle, and high tertiles (left, center, 
and right columns). The emission lines that are fit in these stacked spectra are H/3, the [O III] doublet, Ha and the [S II] doublet. The 
wavelengths of these emission lines are marked with dashed lines. 

the uncertainties in the subsequent fits. Therefore, when 
fitting the emission lines, the interpolated uncertainties 
must be corrected by a factor v^ft, with R representing 
the ratio of the length of the bin in the original spec- 
trum and the length of the bin in the interpolated spec- 
trum. Later, every galaxy is normalized by its Ha flux. 
This step is done to avoid having the stacks be domi- 
nated by the few brightest galaxies. The emission lines 
are closer to the edges of the stacked spectra. Thus, 
different lengths are fit to the left and to the right of 
the emission lines (300 and 600 A) in order to account 
correctly for the small residual continuum. An exam- 
ple of the stacked emission-line fits is shown in Figure [3j 
Emission-line fluxes are then extracted from these com- 
posite spectra (follo wing the emission- line extraction pro- 
cedure described in j ]2.3[ ), which will represent statistical 
properties of the galaxies in every Ha luminosity, stellar 
mass, and Ha equivalent width bin. All the uncertainties 
are given by the p ython version of the mpf it fitting code 
( Markwardt 2009 ) but we checked that those uncertainty 
values are compatible with uncertainties derived from a 
Monte Carlo simulation. 





(mOII 



Tiiiii) 



Figure 3. An example of the fits to the stacked spectra. The 
emission-line fits are plotted with red lines. The wavelengths of 
the emission lines are marked with dashed lines and the observed 
Lhq bin is annotated in a box. 

The number of data points in the re-sampled spec- 
tra needs to be considered in order to correctly estimate 
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Figure 4. The observed Ha luminosity versus stellar mass of our 
galaxy sample. 

There is the concern that the Balmer decrement de- 
termined from stacked spectra may not give the average 
extinction in that luminosity bin. The stacked result is 
dependent upon the intrinsic distribution of individual 
Balmer decrements, whether it is flat, gaussian or, more 
likely, gaussian with a tail to higher extinction. To test 
this potential bias, we ran a Monte Carlo simulation to 
make stacks from individual SDSS galaxies culled from 
the larger Balmer decrement distribution (in a given lu- 
minosity bin). We then measured the Balmer decrement 
in the stack and used the inferred extinction to determine 
the total intrinsic Ha luminosity density from all of the 
galaxies used in the stack. The total Ha flux was com- 
pared to the total flux when summing the Ha fluxes cor- 
rected for extinction based on individual Balmer decre- 
ment measurements for each galaxy. The results indicate 
that these two methods agree within ±20%, despite the 
fact that the distribution of Balmer decrements changes 
significantly from low to high Lhq in SDSS. These uncer- 
tainties are within the uncertainties of the Balmer decre- 
ment analyzed in the present work, so we do not consider 
them a significant bias. 

2.5. Galaxy stellar mass from population synthesis 
models 

Galaxy stellar mass is the most robust physical pa- 
rameter derived from spectral-energy-distribution (SED) 
fitting. This technique has been widely used in the lit 
erature (e.g., [Sawicki fc Yee||1998[ |Papovich et al. 2001 
Shapley et al.|2i 



2( )01| |Bundy et al.l2O06|]Barro et al |201l| ) 



In this section, we calculate stellar masses tor all galax- 
ics in our sample fitting multiwavelength photometry 
(F475X, F600L P, F110W, F160W and IR AC 3.6 ^m; see 
Spl) to a grid of Bruzual & Chariot ( 2003 ) st ellar p opula- 
tion synthesis models assuming a|CEabrier ( 2003 ) initial 

: et al.||20 



mass function. The FAST code ( |Kriek et al.||2009| ) is 
used to derive these stellar masses! The FAST input 
parameters are: the stellar optical extinction Ay cnaT is 
set to the range — 4 mag in steps of 0.1 mag and the 
metallicity by mass fraction is fixed to 0.008 (solar metal- 
licity is 0.02 in these units). The age of the galaxies is 
constrained in the range between 10 7 — 10 10 yr and the 
star-formation history is modeled as an exponentially de- 
clining tau model exp(— t/r) with log 10 (r/yr) = 7 — 10 
in steps of 0.1. The contribution of the emission lines to 



the broad-band flu xes is removed b efore the SED fitting 
(as we suggested in Atek et al.|2011 ) . The Ha luminosity 
(uncorrected for dust absorption) versus stellar mass for 
galaxies in our sample is shown in Figure |4j where more 
luminous Ha emission is seen in more massive galaxies. 
Figur e [4] is related to the star- format ion main sequence 



(e.g., Noeske et al. 2007 Dommguez et al., in prepara- 
tion). 

2.6. Contamination of Balmer line measurements 

Our measurement of the Balmer emission line fluxes 
may be significantly affected by strong stellar absorption 
lines, unresolved emission lines near the Balmer lines, 
and contributions from AGNs. Below we estimate the 
magnitude of each of these effects. 

2.6.1. Stellar absorption lines 

Absorption lines from stars in the galaxy may affect 
the emission-line measurements. We estimate the Ha 
and H/3 absorption lines from the best stellar popula- 



tion model from the SED fitting described above \ 



2.5 



The Ha and H/3 equivalent width (EW) distributions 
are shown in Figure [5] for the three galaxy mass bins 
that were shown in the middle panel of Figure [2] From 
the lowest to the highest mass bin, the median and stan- 
dard deviation of the Ha EWs in absorption distribu- 
tions are (2.42±0.51) A, (2.72±0.52) A, (3.05±0.67) A. 
The corresponding values of H/3 EWs in absorption are 
(3.22±0.82) A, (3.60±0.87) A, (4.09±1.15) A. We com- 
pute the median of the distribution of the individual Ha 
and H/3 EWs in emission that contribute to each stack. 
Then, the ratios between these medians and the EWs in 
absorption corrections are estimated. Since the fluxes are 
proportionally related with the EWs, we use these ratios 
to convert to flux corrections independently both for Ha 
and H/3. The uncertainties in the absorption corrections 
are propagated throughout our analysis. 

2.6.2. The [N II] contribution to Ha 

The resolution of the WISP spectra is not high enough 
to resolve the [N II] A 6548, 6583 doublet from the Ha 
line. This may bias our estimate of the Ha flux to- 
wards higher values, which may lead to systematically 
high Balmer decrements. Of particular concern is that 
the [N II] contribut ion to Ha is expe cted to increase with 



stellar mass (e.g., Erb et al. 2006). Furthermore, the 



Ha/[N II] ratio at a given mass evolves with redshift 



HErb et al-pPOu] |Sobral et~a^[2012 ) 

in order to correct tor this effect, we us e the average 



of the expected Ha /JN II] observed locally (Sobral et al 
2012|) and at z ~ 2 i jErb et al.||2006| ). First, we estimate 
the [N II] to Ha contribution at z ~ 2 — 3 by interpo- 
lating between the Ha/[N II] ratios for differen t stellar 
masses provided in Table 2 of Erb et al. ( 2006 ) . These 
fluxes were derived for star-forming g alaxies with simi- 
lar physical properties as our galaxies. Erb et al. ( 2006 1 
found that the [N II] contribution relative to Ha increases 
with stellar masses. For each of our stacks, the geomet- 
ric mean of the stellar masses of the stacked galaxies is 
calculated. Then, from the previous interpolation a cor- 
rection for our observed Ha fluxes is derived. The lErbl 
et al.|(|2006[) ratios consider only the [N II]A6583 line. We 



add an extra correction to also account for the fainter 
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Figure 5. The Ha and H/3 absorption- line distributions of the gaS 
in the text. 

(by a factor of 3) [N II]A6548 line. The second step is 
to estimate the higher [ N III to Ha contrib ution from 



2345 6 2345 6 

Rest-frame equivalent width [A] Rest-frame equivalent width [A] 

axies in our sample estimated in three bins of stellar mass as described 



the relation provided by Sobral et al 
with local star-forming galaxies 



([2012]) calibrated 
i'his relation estimates 



Ha/[N II] as a function of Ha EW. We proceed estimat- 
ing the geometric mean of Ha EW in our stacked spectra 
and apply the local correction. Finally, the total correc- 
tion that we apply for the [N II] contamination of Ha is 
calculated as the average of these two corrections, that is 
from the Erb et al. (2006, approximately 0%, 8%, 19% 
for each stellar mass bin from the least to the most mas- 
sive bin) and Sobral et al. (2012, 9%, 17%, 20%). The 
combined [N II] correction has little effect on the ulti- 
mate results in this work. For instance, for the stacks 
binned in stellar masses these corrections are 4%, 12%, 
19% from the least to the most massive bin respectively. 

2.6.3. Active galactic nuclei contamination 

We investigate whether a significant contribution from 
active galactic nuclei (AGN) may bias our Balmer decre- 
ment calculations. Figure [6] shows [O III]/H/3 versus 
[S II] /Ha for our stacked spect ra. This figure is the BPT 
diagram (Baldwin et al. 1981[ ). The black line s eparates 
the star-form ing and AGIN region as modelled by |Kewley| 
et al. (2001 1. We see in Figure[6]that our stacks are com- 
patible witn the low-metallicify tail of the star-forming 
region rather than the AGN region from data in the local 
Universe from the SDSS DR&p] ( jEisenstein et al.||201l| ). 
Figure [6] also shows for comparison the data from two 
lensed star- forming galaxies located at z ~ 2 (Hainline 
al.||2009[). In suppo rt of these results, the analysis 
py |Sobral etaT] ( [2012] ) shows around 8% of AGN con- 
Ha selec 



let al.||2009t 
by | Sobral e 

tamination in Ha selected samples and higher in more 
luminous galaxies. Since the galaxies presented in the 
present work are typically fainter, the AGN contamina- 
tion is probably even smaller. These results indicate that 
AGN contamination will not significantly affect the re- 
sults of our analysis. The emission-line ratios shown in 
Figure [6] are listed in Table [T] 

3. THEORETICAL BACKGROUND 

The theoretical framework utilized in the subsequent 
sections is described here. Following the empirical ex- 

11 http://www.sdss3.org/dr8/ 
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Figure 6. BPT diagram that distinguishes between star-forming 
galaxies and galaxies with AG Ns. The black line separates the 
star-forming and AGN regions | |Kewley et al.|2001[ |. Our 9 stacks 
are shown with crosses. The cross size is correlated with the bin of 
Ho luminosity, stellar mass, and H« equivalent width (the larger 
the marker the larger the binning quantity). For compariso n, two 
confirmed l ensed star- forming galaxies detected at z ~ 2 by |Hain^] 
|line et al.| j2009| | are shown with a cyan square and a magenta 
circle, Tne contours are from SDSS data. 

tinction relation found in Calzetti et al. (199^^1 all 
quantities in what follows correspond to nebular emission 
lines, not continuum) the intrinsic luminosities, Li ntl are 
given by, 

L mt {\) = L obs (A)10 04 ^ 

= L obs {m nAk(X)E{B - v) (i) 

where L obs are the observed (extinguished) luminosities, 
A\ is the extinction at the wavelength A, and k{\) the 
reddening curve. In what follows, we will use the redden- 
ing curve fc(A) found in COO for our analysis. The color 
excess E(B — V) is defined by, 



E{B -V) = {B — V) obs -(B- V)i 



(2) 



which is the change in the (B — V) color due to dust 
attenuation (that is, the difference between observed and 



12 See also |Witt fc Gordon| pOOo) ; |BeIl fc Kennicutt] \200l\ for 
their concerns about tms empirical approacn. 
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Table 1 

Emission-line ratios of stacked spectra in 0.75 < z < 1.5 
<log 10 X) ([Sll]A6717 + A6732)/HaA6563 [O III]A5007/H/3A4861 



Stacks binned in Ha luminosity, X = Lhq/ ergs" 



-l 



41 57+ - 48 
11.01 _ Q 23 

4l.9lt° ; 1 13 
42.33l L § 9 5 



0.14 ±0.04 
0.14 ± 0.03 
0.23 ±0.02 



4.65 ± 1.10 
2.55 ±0.66 
3.37 ± 1.18 



Stacks binned in galaxy stellar mass, X = M*/Mq 



^+1.38 
3 -0.58 
,+0.35 
-0.34 
10 46+ - 61 



8.56Z 
9. 51"" 



0.08 ±0.03 
0.28 ±0.03 
0.21 ± 0.03 



5.58 ± 1.32 
1.95 ±0.44 
2.15 ±0.92 



Stacks binned in Ha equivalent width, X = EWh^/A 



1.73 
2.07 
2.47 



+0.51 
0.22 
+ 0.11 
0.11 
+0.28 
0.86 



0.24 ± 0.03 
0.26 ±0.03 
0.04 ±0.03 



2.58 ±0.84 
2.91 ± 0.69 
5.89 ± 1.63 



Note. 
Figure [6] 



These are the values for plotting the BPT diagram shown in 



expected color in the absence of dust). The nebular color 
excess is related to the stellar color excess as discussed 
in e.g., Calzetti et al. (2000, hereafter COO). 

We are ultimately interested in deriving dust extinc- 
tions from Balmer decrements. The relationship between 
the nebular emission line color excess and the Balmer 



decrement is give n by (see e.g., the appendix in Mom- 
cheva et al.||2012 for details): 



E(B - V) 



2.5 



E(Hf3 - Ha) 
k(X H p) - k(X Ha ) 



k(X H fi) - k(X Ha ) 



log 



10 



(Ha/H(3) obs 
(Ha/Hf3) mt 



(3) 



where k(Xnp) and fc(Ajj a ) are the reddening curves eval- 
uated at H/3 and Ha wavelengths, respectively. The fac- 
tor E{Hj3 — Ha) is analogous to the color excess but 
defined for H/3 and Ha, instead of the B and U-bands. 
Then, (Ha/ 'H (3) b s is the observed Balmer decrement 
and (Ha/ H f3) int is the intrinsic or unreddened Balmer 
decrement, which is calculated theoretically. The in- 
trinsic Balmer decrement remains roughly constant for 
typical gas con ditions in star-forming galaxies (see Os 



terbrock 1989). In our analysis, we assume the value 



of (Ha/ H /j) int — 2.86, corresponding to a temperature 
T = 10 4 K and an elec tron density n e = 10 2 cm' 3 for 



Case B recombination (Osterbrock 1989). This choice 
is standard for star- forming galaxies m the literature. 
Thus, the nebular color excess is given by the following 
equation, 



E(B-V) = 1.97 log 



10 



(Ha/HP) 



obs 



2.86 



(4) 



Coming back to Equation [T] A\ is the extinction in 
magnitudes at the wavelength A, 



A\ = k(X)E(B — V). 



(5) 



In particular, selecting the COO reddening curve, Equa- 
tion [5] turns into, 



A Ha = (3.33 ± 0.80) x E(B - V) 



A v = (4.05 ± 0.80) x E(B - V) 



(6) 



(7) 



which are the Ha extinction and the optical extinction, 
respectively. The parameters E(B — V), An a , and Ay 
will be given in Table [2] for our stacked spectra. Any po- 
tential effect coming from galaxy inclinations should be 
alleviated from our stacking methodology since we are 
averaging a statistical sample of galaxies with indepen- 
dent inclinations. 



4. RESULTS 

4.1. Balmer decrements from Ha luminosity stacks 

Balmer decrements as a function of Ha luminosities 
(Z/Ha) a re analyzed in this section. The stacked spectra 
in three different observed Ln a bins are shown in the top 
panel of Figure [2j The total number of stacked galaxies 
is 128. The size of bins are chosen to have the same 
number of galaxies. 

Figure [7] (left panel) shows the Balmer decrement as 
a fun ction of observed Ha luminosity. As discussed in 
{2.6 the absorption- line corrections were applied to the 



Ha and H/3 fluxes and also a correction to account for 
the contribution of the [N II] doublet to Ha. These 
corrected Balmer decrements are shown with filled red 
stars. The open triangles show the Balmer decrements 
before applying any of these corrections and the open 
squares show the Balmer decrements after applying only 
the absorption-line corrections. The uncertainties are 
omitted in these open symbols for clarity. The intrinsic 
Balmer decrement considered in this work of 2.86 (dust- 
free galaxies) is shown with a dotted line. For compari- 
son, star-forming galaxies in the local Universe from the 
SDSS are shown with contours (galaxies with AGNs were 
not removed from this sample). The relation between 
Balmer decrements and observed Lna for these galaxies 
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Figure 7. Balmer decrements versus observed Ha luminosity (left panel), galaxy stellar mass (middle panel), and rest-frame Ho equivalent 
width (right panel). In all the panels, the open red triangles show the stacked spectra from the whole redshift range without correcting 
for Ha and H/3 line absorption or [N II] contamination. The open red squares are the Balmer decrements that include the absorption-line 
corrections (the uncertainties of the open symbols are omitted for clarity). The filled red stars show the stacked spectra from the w hole 
redshift range 0.75 < z < 1.5, which are corrected for both Ha and H/3_ absorption lines and the [N II] doublet contamination (see j ]2.6| 
for details). All the values represented in this figure are listed in Table ^1 The dashed-magenta line in the left pan el shows the results by 
|Hopkins et al.|l|2001[ |. A dashe d-blue line in t he mi ddle panel shows the relationship provided by Garn & Best (2010) for local galaxies. The 
results at z ~ 1.5 presented by |Sobral et al.| l |2012| l are shown with green circles in the left and middle panels. Uontours show star-forming 
and AGN SDSS galaxies for comparison. The dotted-black line shows the intrinsic Balmer decrement assumed in our analysis of 2.86. 
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Figure 8. Balmer decrement versus Ha luminosity for individual 
spectra. The black crosses represent galaxies with Ha/H/3 detected 
at more than 3<r (11 galaxies). The gray arrows show 3<r lower 
limits. Other symbols are the same as in the left panel of Figure[7] 



from the analysis presented in Hopkins et al.| ( 2001 ) is 
shown with dashed-ma genta line. The SPSS contours do 



not exactly match the Hopkins et al. (2001) results be- 



ca use the samples a re not exactly the same. The results 



by |Ly et ah ([2012 1 at z ~ 0.5 are not explicitly shown 
in Figure jj] (left panel) since they do not pr ovide any 
param eterization. At any rate, the results by |Ly et aL 
( |2012[ ) are compatible within uncertainties with the re- 
sults presented in |Hopkins et al. ( 2001 ). The green circles 
in Figure [7] are the results found in Sobral et al. (2012) 
at z ~ 1.5~Table[2]lists the main dust properties of the 
stacked spectra. 

Figure M shows the distribution of the individual galax- 
ies m Ho/H/3 versus Ha luminosity. The black crosses 
show galaxies with Ha/H/3 detected at more than 3a 
whereas the gray arrows represent 3a lower limits. There 
is a broad scatter in this graph for a given observed Ha 
luminosity that suggest that our extinction results can 
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Figure 9. The [O III] /Ha ratio versus Ha luminosity. The dust- 
corrected data have been slightly shifted in the a;-axis for clarity. 
The dashed line represents the constant value of unity. 

only be applied to statistical galaxy samples. 

For the faintest star-forming galaxies [O III]A5007 is 
brighter than Ha. This fact is seen clearly in Fig- 
ure [9j where we show the [O III]A5007/Ha ratio as a 
function of Ha luminosity. The red squares show the 
observed [O III] A5007/Ha. The blue circles show the 
[O III]A5007/Ha after correcting for dust extinction us- 
ing the Balmer decrements at every observed Ha lumi- 
nosity bin. The trend of increasing [O III]A5007/Ha ratio 
can not be explained by dust extinction since the relative 
corrections between the two wavelengths are generally 
small. After correcting for dust extinction, the faintest 
luminosity bin has higher [O III]A5007/Ha > 1, signifi- 
cantly higher than in the higher luminosity bins. The rel- 
ative increase in [O III] A5007 flux is likely due to decreas- 
ing metallicity a nd /or ionization parameter (e.g., |Kewley 
fc Dopita|[2002l ). 



4.2. Balmer decrements from galaxy stellar mass stacks 
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Table 2 

Dust properties for WISP galaxies 





<i°g 10 P0> 


H a /H h H a /Hp c H a /Hp d E(B - V) c [mag] 


A Ha e [mag] 


A v c [mag] 








Galaxy spectra stacked in Ha luminosity, X = Ln^/crgs 


-l 




43 
43 
42 


41 t-y+0.23 

41 Qi+0.13 

42 ^Q+l- 05 


3.99 ±0.95 3.51 ±0.85 3.29 ± 0.82 0.12 ±0.21 
4.14 ± 1.03 3.59 ±0.91 3.08 ± 0.79 0.06 ± 0.22 
8.84 ±2.66 5.84 ± 1.86 5.01 ± 1.60 0.48 ± 0.27 


0.40 ±0.71 
0.21 ± 0.73 
1.59 ±0.98 


0.48 ±0.86 
0.25 ±0.89 
1.94 ± 1.17 


Galaxy spectra stacked in stellar mass, X = M* /M© 


43 
43 
42 


8.56±°;gj 
10.46lJ;« 


3.92 ±0.93 3.66 ±0.87 3.49 ± 0.85 0.17 ±0.21 
5.16 ± 1.10 3.99 ±0.88 3.50 ± 0.78 0.17 ±0.19 
6.80 ± 2.79 5.52 ±2.29 4.46 ± 1.86 0.38 ± 0.36 


0.57 ± 0.71 
0.57 ± 0.65 
1.26 ± 1.22 


0.69 ± 0.85 
0.70 ± 0.78 
1.54 ± 1.47 




Galaxy spectra stacked in rest-frame Ha equivalent width, X = 


EWhc/A 




43 
43 
42 


7 o+0.22 
'-■' -Q.Bl 
2 07+ ' 11 
9 A7 +0.86 
z -*' -0.28 


6.93 ±2.16 3.86 ± 1.36 3.18 ± 1.13 0.09 ± 0.30 
3.96 ±0.90 3.46 ±0.80 3.05 ± 0.72 0.06 ± 0.20 
5.69 ± 1.58 4.88 ± 1.37 4.72 ± 1.35 0.43 ± 0.24 


0.30 ± 1.01 
0.19 ±0.67 
1.42 ±0.88 


0.36 ± 1.23 
0.23 ±0.81 
1.73 ± 1.05 


Note. 


— All observables 


are given for nebular properties. 







a N is the number of stacked galaxy spectra in the bin 
b No correction applied. 



c Corrected only for Ha and H/3 absorption lines (see |2 

Corrected for Ha and H/3 absorption lines and [N II] contamination (see f 2.6 i 
c Calculated from the absorption line and [N II] corrected Balmer decrements. 



We also analyze the Balmer decrement from the 
stacked spectra versus galaxy stellar mass. Figure [7] 
(middle panel) shows the Balmer decrements calculated 
for the full galaxy sample. (The three composite spectra 
are shown in the middle panel of Figure^ as well.) The 
symbols are the same as in the left panel of Figure [7} 
We also show star-forming galaxies in the local Universe 
from the SDSS with contours (as in the previous section, 
galaxies with AGNs were not removed) . The dashed-blue 
line shows the relationship between the Balmer decre- 
ment and galaxy stellar mass calc ulated from the anal - 
ysis of a different SDSS sample by |Garn fc Best] fl2010|). 
The g reen circles are the results found in Sobral et al. 



(2012) at z ~ 1.5. Table[2]lists the dust properties of the 
galaxy spectra stacked in stellar masses. 

4.3. Balmer decrements from Ha equivalent width 
stacks 

The equivalent width of the Ha emission line (EWh q ) 
is a proxy for the specific SFR of a galaxy. This quan- 
tity is defined as the width of the integrated continuum 
underneath the line that matches the flux of the line. 
The rest-frame Ha EW are used in this section, this is 
the observed Ha EW divided by (1 + z). We study as 
well the Balmer decrement as a function of EWh q - The 
stacked spectra for the overall sample, binned in EWh , 
are shown in the bottom panel of Figure [2] The right 
panel of Figure [71 shows the Balmer decrement calcu- 
lated from the WISP stacked spectra as a function of Ha 
EW. All colors and symbols are the same as in the other 
panels of Figure fjj The main properties of the Ha EW 
based stacks are listed in Table [2j 

5. DISCUSSION 
5.1. Balmer decrements at 0.75 < z < 1.5 



We find from Figure [jj (left panel) that the Balmer 
decrement, and therefore dust extinction, increases with 
observed Lh q in the luminosity range of our analysis. 
The increase of the Balmer decrement with observed L-^ a 
was a l ready e stablished at lower r e dshifts by | Brinchmann| 
etldT1 ( [2004| ; |Garn fc Be^ pOTOj ) ; |Ly et al.|l20T2 [) This 
trend was also found recently by Sobral et al.| (20121 
at z ~ 1.5 using le ss robust dust indic ators than the 
Balmer decrement. Sobral et al. ( 2012[ ) used the ratio 
Ha/[0 II] (taken with narrow-band photometry), which 
is dependent on the unknown galaxy metallicity. We also 
note that the high star-forming galaxies represented in 
our sample are found rarely in the local Universe. On 
the other hand, galaxies with the lowest observed Lh q 
(10 41 erg s~ 4 < Lhq ^ 10 42 erg s _1 ) were found in §4]to 
have Balmer decrements compatible with dust-free galax- 
ies. 

Our data are not enough to characterize the depen- 
dence of the Balmer decrement with galaxy stellar mass 
but it is consistent with no evolution of the local relation 
(Figure [71 middle panel). A trend of the Balmer decre- 
ment with EWhq (which traces specific star formation) 
is also not evident from the right panel of Figure [7J 

5.2. Evolution of the Balmer decrement with redshift 

There is a redshift evolution in dust properties as a 
function of observed Lhq from the local Universe to 
higher redshifts. This result remains valid even without 
applying the [N II] corr ections. Our results are compat- 
ible with the claims by Sobral et al. ( 2012[ ) that an 
galax^plat z ~ 1 (between 10 4 ^° andTF 250 erg s"\ 



Sobrafet al. 

ga laxy today (approximately 10 4120 



2007) 



2012|) has the same dust extincti on as an 

erg s _1 , Ly et al. 

This result comes from the fact that the evolution 



13 According to a Schechter 1 19761 parameterization 
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in the dust-extinction dependence with Lhq is similar to 
the evolution from the local Universe to z ~ 1 of X* in 
the luminosity function. On the other hand, looking at 
the middle panel of Figure [7J we see that our galaxies 
span the same stellar mass range as the galaxies studied 
in the local Universe by the SDSS. The dependence of the 
Balmer decrement with stellar masses does not seem to 
evolve with time, w hich is in agreement with results by 
Sobral et al. (2012) utilizing photometry and a different 
tracer of extinct ion . 

These results are consistent with the following picture 
of galaxy evolution. Young star-forming galaxies have 
typically low content of dust and metals. Significantly 
older galaxies that built many stars in the past will be 
dustier and metal rich. When observing younger galaxies 
at earlier cosmic times, we see less dust in the interstel- 
lar m edium even at fixed SFR (see e.g., Reddy et al. 
20i0| 



Universe up to 



1.5 as a function of Lh q is found 



Burgarella et al. 


20071) or less ro 


dicators (e.g., 


So 


bral et al.|2012 



2006 



Nevertheless, the galaxy stellar mass is cor 
related with the production of dust and metals in such 
a way that the same amount of dust is produced per 
unit of formed stellar mass regardless of which redshift 
it was formed. Typically the dust is mainly produced by 
asymptotic giant branch (AGB) stars, which are formed 
when the galaxy is several 100 Myr old. Therefore, we 
expect that the correlation of dust extinction with stel- 
lar mass will be similar at any epoch except at higher 
redshifts, when most stars have not had time to enter 
the AGB phase. Not surprisingly, we find that this issue 
is not significant at the epoch around z ~ 1, when the 
universe was nearly half its present age. However, there 
is one caveat (which is still controversial) that we do not 
consider in this picture: dust may be destroyed and/or 
expelled from the galaxy by outflows. 

A possible redshift evolution in the Balmer decrements 
within our sample is investigated by means of binning in 
two different redshifts ranges but no conclusion can be 
drawn from our sample due to large uncertainties. 

We note that estimating dust extinction from the 
Balmer decrement is based on the assumption that the 
Ha extinction is constant across all Ha emitting region. 
As with any measure of extinction derived from unextin- 
guished light, the Balmer decrement can not account for 
optically thick regions of the galaxies. This ma y result 
in an unde restimate of the total extinction (e.g., Meurer 
et al. 120021). However, because the extinction is higher in 



the ultraviolet and because the intrinsic ratios are known, 
we believe that our Balmer line-derived estimates of the 
extinction are more robust than those derived with ultra- 
violet colors. Also, when comparing to other epochs to 
determine an evolution in extinction, we use the same ex- 
tinction metrics for consistency (the Balmer decrement). 
Ideally, we would like to also measure the infrared lu- 
minosities of these galaxies to fully account for obscured 
star formation, but that is not possible given the resolu- 
tion and confusion limits of existing infrared telescopes. 

6. SUMMARY 

Statistical interstellar-dust properties of star-forming 
galaxies were analyzed in this work as a function of Ln a , 
galaxy stellar mass, and rest-frame Ha equivalent width. 
These properties were derived by stacking the spectro- 
scopic observations of Balmer decrements of 128 star- 
forming galaxies from the WISP survey. 

Evolution in dust-extinction properties from the local 



more attenuated in the local Universe than at z ~ 1 by 
a factor that is dependent on luminosity. However, the 
Balmer decrement as a function of galaxy stellar mass 
does not seem to evolve with redshift (as suggested by 
Sobral et aO2012| ). 

We show lhat the fainter star-forming galaxies have 
brighter [O III] than Ha emission lines. Future galaxy 
surveys of faint star-forming galaxies may take advantage 
of the fact that [O III] is typically brighter than Ha at 
Lua < 5 x 10 41 erg s _1 as it will allow for detection of 
fainter galaxies and/or higher redshifts than Ha-targeted 
surveys. 

The number of sources in the WISP catalog will in- 
crease considerably before its completion and will benefit 
from ongoing ground-based observations. Our dust ex- 
tinction results will be updated when these observations 
are available. The WISP team will soon present further 
results on the Ha luminosity function (Colbert et al., in 
preparation), massive galaxies (Bedregal et al., in prepa- 
ration), mass-metallicity relation (Henry et al., in prepa- 
ration) and star-forming main sequence (Dommguez et 
al., in preparation). These results will allow us to present 
a comprehensive picture of the Universe at z ~ 1 and be- 
yond. 

The authors thank Naveen Reddy for reading and im- 
proving the paper and David Sobral for fruitful discus- 
sions. This work is based on observations made with 
the NASA/ESA Hubble Space Telescope, which is op- 
erated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with programs 11696, 
12283,12568. 
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